function y = rst_sta(temp, rym)
    Fs = 8192;
    wavef = [523, 587, 659, 698, 783, 880, 988];
    x = linspace(0, 2 * pi * rym, floor(Fs * rym));
    y = sin(wavef(temp) * x);
end

